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ABSTRACT 

A mathematical framework is presented for use in the experimental determination of the polarimetric 
response of observatory instrumentation. Elementary principles of linear algebra are applied to model 
the full matrix description of the polarization measurement equation by least-squares estimation of 
non-linear, scalar parameters. The formalism is applied to calibrate the center element of the Parkes 
Multibeam receiver using observations of the millisecond pulsar, PSR J0437— 4715, and the radio galaxy, 
3C 218 (Hydra A). 

Subject headings: methods: data analysis — instrumentation: polarimeters — polarization — pulsars: 
individual (PSR J0437-4715) — techniques: polarimetric 
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1. INTRODUCTION 

Polarization measurements provide additional insight 
into the phenomena involved in both the emission and 
propagation of electromagnetic radiation. However, the 
processes of reception and detection introduce instrumen- 
tal artifacts that must be corrected before meaningful in- 
terpretations of experimental data can be made. Ideally, 
the instrumental response is estimated by observing at 
least two calibrator sources of well-determined, polarized 
radiation (Hamaker, Bregman & Sault 1996). Under the 
assumption that the observatory apparatus respond lin- 
early and remain stable between calibrator observations, 
the known instrumental response is then inverted and used 
to calibrate observations of other sources. 

In the absence of a sufficient number of calibrator 
sources, instrumental calibration may in some cases be 
performed by fitting the available polarimetric data to a 
predictive model. This model must describe the observed 
sources of radiation, the response of the instrument, and 
any additional propagation effects that arise in the inter- 
vening media through which the signals are transmitted, 
such as the Earth's ionosphere. In order for the model to 
be constrained, at least one component of the applied mea- 
surement equation must vary as some function of one or 
more independent variables. For example, the constrain- 
ing transformation may consist of the geometric projection 
of the antenna receptors onto the plane of the sky, which 
may vary during transit of the source for antennas without 
an equatorial mount. 

For an altitude-azimuth mounted antenna, the receiver 
feed is rotated about the line of sight through the paral- 
lactic angle, a constraint that has been utilized in a num- 
ber of previous studies, including Stinebring et al. (1984), 
Turlo et al. (1985), Xilouris (1991), McKinnon (1992), and 
Johnston (2002). In each of these matrix prod- 

uct is expanded to yield a set of scalar equations that 
describe the measured Stokes parameters as a function of 
parallactic angle; these are then solved using conventional 
techniques. This approach has two significant limitations. 
First, in order to simplify the derivation of the scalar equa- 
tions and their partial derivatives, it is necessary to make 
small value approximations that are not generally valid. 



Second, as any alteration of the reception model neces- 
sitates a laborious expansion of the matrix product, this 
approach is non-conducive to experimentation with a va- 
riety of parameterizations. 

In addition to these shortcomings, other fundamental 
limitations and conceptual errors have been incorporated 
into previous treatments. For instance, Turlo et al. (1985) 
omit circular polarization from consideration, an over- 
simplification that restricts further application of their re- 
sult to antenna with perfect, circularly polarized recep- 
tors. Stinebring et al. (1984), followed by Xilouris (1991) 
and McKinnon (1992), begin with the assumption that 
the complex gains of the two polarizations may be inde- 
pendently calibrated prior to determination of the cross- 
coupled antenna response. However, this approximation is 
accurate only when the differential gain and phase trans- 
formations commute with the antenna response matrix, 
a condition satisfied only when the feed receptors have 
well-determined, orthogonal polarizations. In general, the 
complex gains may be accurately quantified only after the 
effect of the antenna feed is included. 

This error is essentially compounded in Johnston (2002) 
where, although the instrumental response parameters are 
jointly estimated, the measured data are first corrected 
using separately-determined, inaccurate estimates of dif- 
ferential gain and phase. In addition, the observed Stokes 
parameters are incorrectly normalized by the total inten- 
sity; a quantity such as the invariant interval must be used 
(Stinebring et al. 1984; Britton 2000). Finally, as only ob- 
servations of an unknown source are used to constrain the 
instrumental response parameters, the solution derived by 
Johnston (2002) is not unique. Rather, as shown in Ap- 
pendix B, there exists a degenerate set of solutions that 
can be resolved only with additional observations of one 
or more sources with known circular polarization and po- 
sition angle. 

In contrast to previous work, the treatment presented in 
this paper is based entirely on a full matrix description of 
the polarization measurement equation. All of the matrix 
products, including those of the required partial deriva- 
tives, are evaluated in software; therefore, no small value 
approximations are necessary and a variety of model spec- 
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ifications may be tested with relative ease. Furthermore, 
all of the model parameters are jointly constrained using 
observations of both unknown and partially known source 
polarizations. Following a brief mathematical review in 
Section 2, the method employed to solve the matrix equa- 
tions is developed in Section 3. Two possible decompo- 
sitions of the instrumental response are described in Sec- 
tion 4, one of which is applied in Section 5 to the calibra- 
tion of the Multibeam receiver and downconversion system 
utilized by CPSR-II: the 128 MHz baseband recording and 
real-time processing system at the Parkes Observatory. 

2. POLARIZED RADIATION AND PROPAGATION 

Electromagnetic radiation is described by the two 
complex-valued components of the transverse electric field 
vector, e(t) = (eo(t), e\ (t)). The measurable properties 
of e(t) are represented by the coherency matrix, p = 
(e(i)® e^(t)), where the angular brackets denote time av- 
eraging and is the Hermitian transpose of e. The co- 
herency matrix may be written as a linear combination of 
Hermitian basis matrices (Britton 2000), 



{S <t + S • <x)/2, 



(1) 



k=0 



where <x is the 2x2 identity matrix, er = (eri, <r 2 , 0-3) 
arc the Pauli spin matrices, So is the total intensity, or 
Stokes I, and S = (Si, S2, S3) is the Stokes polarization 
vector. As described in Appendix A, S = (Q, U, V) in 
the Cartesian reference frame applied in this paper. The 
Stokes parameters may also be expressed in terms of the 
coherency matrix (Hamakcr 2000), 

S k = tr(<T fc p), (2) 

where tr(A) is the trace of the matrix, A. 

The propagation and reception of the electric field is 
represented by the transformations of linear time- invariant 
systems. Presented with the input signal, e(t), the output 
of a system is given by the convolution, e'(t) = j(t) * e(t), 
where j(i), is the complex- valued, 2x2 impulse response 
matrix. By the convolution theorem, this transformation 
is equivalent to E'(w) = 3(oj)E(u>), where J is the familiar 
Jones matrix. When observed over a sufficiently narrow 
band, the variation of J with frequency is ignored, and po- 
larimetric transformations are represented by Jones matrix 
multiplications in the time domain. Under the operation, 
e'(t) = Je(t), the coherency matrix is transformed as 



p' = JpJ f . 



(3) 



This congruence transformation is called the polarization 
measurement equation (Hamaker 2000). It forms the basis 
through which measured quantities are related to the in- 
trinsic polarizations of the sources and used to model the 
unknown instrumental response. 

3. MAXIMUM LIKELIHOOD ESTIMATOR 

To solve the polarization measurement equation using 
conventional methods of least-squares minimization, it is 
necessary to design a scalar figure-of-merit function and 
to calculate both its gradient and curvature with respect 
to scalar model parameters. Let r] represent the vector of 
scalar parameters that describe the model, including the 



instrumental response and the polarizations of the sources. 
Furthermore, let x represent the vector of independent 
variables that constrain the measurement equation, such 
as the observing frequency and epoch. Given x, the model, 

f/ m (x;r,); l<m<M, (4) 

must predict the measured polarization of each source, 
where M is the number of sources. In the interest of 
brevity, the model may also be represented by p' m (x). 

Now consider N m independent observations of the m th 
source, each made at a unique coordinate, x ni ni to yield 
the measured Stokes parameters and their estimated er- 
rors, {S' k ± <Jk}m,n- The best-fit model parameters will 
minimize the objective merit function, 

M N m 3 ,™ _ ci / ni 2 

2/ \ V* \-°k,m,n J k,m\^m,m '1)1 , . 

x(v) = 2^1^1^^^ — — ( 5 ) 

m— 1 n— 1 k— 



2 

^k,m,n 



where S' ktm (x m!n ;rj) = tr[er fe p' m (x m . n ; -q)] are the Stokes 

parameters of the m th source as predicted by the model 
at x mtn . The gradient of x 2 with respect to the scalar 
parameters, rj, has components, 



where 



.(»?) 



k=0 



(7) 



k.m.n 



(8) 



Taking an additional partial derivative yields 

d 2 X 2 = _ 2 V Vtrf dA m ,„ dp' m (x m . n ) 
dri s dr] r V d-q s dr\ r 

m—l n—l v 

where, following the discussion in § 15.5 of Numerical 
Recipes (Press et al. 1992), the term containing a sec- 
ond derivative in equation (8) has been eliminated. Us- 
ing equations (5) through (8), the Levenberg-Marquardt 
method is applied to find the parameters that minimize 

x 2 (v)- 

4. PARAMETERIZATION OF THE MODEL 

It remains to specify the scalar values, rj, that param- 
eterize the polarization measurement equation as well as 
the partial derivatives of p' m (x;r]) with respect to those 
parameters. As in equation (1), the coherency matrix, p m , 
of each input source polarization is completely specified by 
the four Stokes parameters, Sk, m - The partial derivatives 
with respect to these parameters are simply 



dS k 



Vk 

2 ' 



Now, let the model of the instrumental response be rep- 
resented by the complex- valued 2x2 Jones matrix, J. If J 
satisfies equation (3) then it belongs to the set of solutions 
given by J(</>) = e^J, where i = \/^T. That is, the co- 
herency matrix is insensitive to the absolute phase of the 
signal and 4> may be arbitrarily chosen, leaving seven de- 
grees of freedom with which to describe the instrumental 
response. Two possible parameterizations are considered: 
the algebraic decomposition employed by Hamakcr (2000) 
and the phenomenological description of Britton (2000). 
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4.1. Algebraic Model 

Following Hamaker (2000), an arbitrary matrix, J, is 
represented by its polar decomposition, 

J = JB A (/3)R ft (0) (10) 
where J = (det J) 1 / 2 and, as described in Appendix A, 

Bm(/3) = cosh [3a- + sinh (3rh • cr, and (11) 
R-n(^) = cos0er o + isin^n • <r. (12) 
The phase of the complex value, J, is set to zero and, as 
any rotation about an arbitrary axis may be decomposed 
into a series of rotations about three perpendicular axes, 
the instrumental response is written as 

3 

J H = GB™(/3) JjR 4h (0 fc ), (13) 

k=l 

where G — \J\ is the absolute gain and Sfc are the or- 
thonormal basis vectors defined in Appendix A. Equa- 
tion (11) may be written as B ? f l (/3) = bo ""o + b • cr, where 
bo = cosh/3 and 6 = sinh /3m. Noting that bo > 1 and 
6q = 1 + 6-6, the three degrees of freedom of Hm((3) 
are specified by 6 = {b\, 6 2 , 63). Therefore, the instrumen- 
tal response is parameterized by G, 61-3, and 4>i-3- The 
partial derivatives of Jh with respect to these seven pa- 
rameters are calculated using products of 

b k 

cr + cr k , (14) 



and 



db k 
dH Sk (4> k ) 



Vl + b-b 



sm<j> k cr + icos<j>k cr k . 



(15) 



4.2. Phenomenological Model 



S a = 



(17) 



Beginning with Britton (2000), the response of an ideal 
feed with two orthogonally polarized receptors is given by 
S(9,e) = Hu(e)Hij(9). Here, the receptors have elliptic- 
ities equal to e and mutually perpendicular orientations 
defined by 9. The basis vectors, q, u, and v, are defined 
in Appendix A. Using this notation, a receiver with non- 
orthogonal receptors is represented by 

C = 5 o S(0o,eo) + *iS(0i,ei), (16) 
where d a is the 2x2 selection matrix, 

S a 

s la 

S a b is the Kronecker delta, and the product, £ a B, returns a 
matrix that contains only the a th row of B. Equation (16) 
is equivalent to equation 16 of Britton (2000) and is used 
without making any first order approximations or further 
assumptions about the feed. 

The differential gain and phase of the instrument are 
represented by 13g 1 (7) and Hg 1 (tp) where, as in equa- 
tions (Al) and (A2), 7 and tp parameterize Lorentz boost 
and rotation transformations (Britton 2000). Including 
the absolute gain, G, the instrumental response is written 
as 

3 B =GB, 1 ( 7 )R, 1 ( V )C. (18) 

The partial derivatives of with respect to its seven 
scalar parameters, G, 7, tp, 9o-i, and eo-i, are calculated 
using equation (15) and 

dB Sh ([3 k ) 



5. APPLICATION 

Radio pulsar observations provide an excellent source 
of data with which to constrain the polarization mea- 
surement equation. The often highly polarized state of 
a pulsar signal can vary significantly as a function of pulse 
longitude. When integrated over a sufficient number of 
spin periods, mean polarimctric pulse profiles generally 
exhibit excellent stability on timescales much longer than 
those over which calibration observations are made. Con- 
sequently, multiple on-pulse longitudes from a single pul- 
sar may be included as unique and stable input source 
polarizations, greatly increasing the number of available 
constraints when compared with non-pulsed sources. Fur- 
thermore, any non-pulsed background polarization is effec- 
tively eliminated by subtracting the off-pulse mean from 
each integrated pulse profile. 

The millisecond pulsar, PSR J0437— 4715, represents an 
ideal candidate for use in the regular calibration of the 
apparatus at the Parkes Observatory. As part of the high- 
precision pulsar timing program, it is often observed from 
rise to set, providing measurements with a wide range in 
parallactic angle. The model that describes these obser- 
vations has the form, 



p^($)=JR e ($)p m R c t(c i >)jt ; 



(20) 



d(3 k 



sinh P k (To + cosh f3 k cr k . 



(19) 



where R# ($) is the rotation about the line of sight by the 
parallactic angle, $. 

As demonstrated in Appendix B, additional calibrator 
observations are required in order to uniquely determine 
the solution to equation (20). The Parkes Multibeam re- 
ceiver is equipped with a noise diode that ideally injects 
a 100% linearly polarized signal with a position angle of 
45 degrees into the receiver feed horn. This reference 
source is switched using a signal generator and the ob- 
served square waveform is integrated modulo its period. 
Ideally, the on-pulse longitudes of the input reference sig- 
nal contain additional flux with Stokes parameters given 
by Go [1, 0, 1, 0], where Go is the reference flux density. Un- 
der this assumption, inclusion of reference signal obser- 
vations breaks the degeneracy by constraining the boost, 
(3 V , which mixes Stokes I and V, and the rotation, (f> v , 
which mixes Stokes Q and U. In a separate flux calibration 
procedure, the reference signal is observed simultaneously 
with the bright Fanaroff-Riley type I radio galaxy, 3C 218 
(Hydra A) , producing absolute flux estimates of both the 
system temperature and reference flux density, Gq. 



5.1. Observations 

Dual-polarization observations of PSR J0437— 4715 and 
Hydra A were made on 19 and 20 July 2003 using the 
center element of the Parkes Multibeam receiver. Two 
64 MHz bands, centered at 1341 and 1405 MHz, were two- 
bit sampled and processed by CPSR-II, the second gener- 
ation of the Caltech-Parkes-Swinburne Recorder. In order 
to maintain optimal linear response during the digitization 
process, the detected power is monitored and the sampling 
thresholds are updated approximately every 30 seconds. 
In addition, the baseband data reduction software corrects 
quantization distortions to the voltage waveform using the 
dynamic level-setting technique (Jenet & Anderson 1998). 



4 



W. van Straten 



Phase-coherent dispersion removal is performed while syn- 
thesizing a 128-channel filterbank (Jcnct ct al. 1997); the 
Stokes parameters are then detected and integrated as a 
function of topocentric pulse phase. Data are averaged for 
five minute intervals, producing uncalibrated mean pulse 
profiles with 2048 phase bins, or an equivalent time reso- 
lution of approximately 2.8 fis. 

The observed flux density varies significantly between 
pulsar and Hydra A observations, resulting in large 
changes to the digitization sampling thresholds. This dif- 
ference is best modeled using the phenomenological de- 
composition of equation (18), applying separate absolute 
gain, differential gain, and differential phase terms to the 
two sets of observations. Also, in order to account for 
phase drifts on timescales of a few hours, the differential 
phase is modeled to vary as a cubic polynomial function of 
time. The three signal path transformations corresponding 
to the three observed sources are summarized in Table 1. 

The flux of the pulsar also varies between observations, 
a result of both intrinsic intensity fluctuations and inter- 
stellar scintillation. Therefore, the measured Stokes pa- 
rameters are normalized by the invariant interval, 

S[ m ^2{detp' m ^ 2 = {S' 2 -S' 2 f' 2 , (21) 

where the polarized intensity, S' = \S |. In order to avoid 
division by small or negative values, those data points with 
S; nv < (To are discarded. Assuming that the measurement 
errors in each of the Stokes parameters are independent of 
each other, the estimated errors in the normalized Stokes 
parameters are given by 

a 2 = S[J S' 2 {S' 2 a 2 + S' 2 a 2 s ), and (22) 

& l = ^inv 4 K^inv + 2S 'k 2 ) a k + <Sfc 2cr i 2 nv], ( 23 ) 

where 1 < k < 3, 

S'^l^S'^al (24) 
fc=i 

and S[ 2 v a 2 m = S^ct 2 + S /2 a 2 . Note that <t is approxi- 
mately proportional to S' when S' <C So- Consequently, 
the normalized total intensity is plotted on a separate scale 
in Figure 1, in which are displayed the normalized Stokes 
parameters from one of the 65 pulse longitudes used to 
constrain the model. 

5.2. Results 

Initial results indicated that the reference signal pro- 
duced by the noise diode of the Parkes Multibeam receiver 
is not actually 100% linearly polarized. Rather, as shown 
in Figure 2, the reference signal consists of ~ 90% linear 
and 3% circular polarization, its position angle is not ex- 
actly 45 degrees (i.e. Stokes Q ^ 0), and it is severely de- 
polarized at the edges of the band due to frequency aliasing 
during downconversion. Consequently, the reference signal 
cannot be trusted as a source with known position angle 
and degree of circular polarization. Also, approximately 
7 MHz from both the top and bottom of the band must be 
discarded as irreversibly corrupted. 

In order to constrain [3 V , it is noted that Hydra A has 
less than 0.1% circular polarization (Roberts et al. 1975). 
Therefore, the off-pulse longitudes of the flux calibration 
observations serve as a source of negligible circular polar- 
ization in the remainder of this development. The un- 
known rotation about the line of sight, cf) v , is artificially 



constrained by assuming that (9o = 0; consequently, all po- 
sition angles are measured with respect to the orientation 
of receptor 0. In modeling the reference signal, Stokes I is 
set to unity, producing an intermediate flux scale in units 
of the reference flux density; Stokes Q, U, and V are varied 
as free model parameters. 

The reception model is solved independently for each 
of the 256 frequency channels in the two 64 MHz bands. 
In each channel are a total of 278 free model parameters, 
corresponding to the Stokes parameters of the 65 selected 
pulse longitudes, the 12 instrumental parameters listed in 
Table 1, and the 6 unknown Stokes parameters shared be- 
tween the noise diode and Hydra A. These parameters are 
constrained by approximately 25220 measured values, pri- 
marily derived from the selected pulse longitudes of the 
97 pulsar observations. Indicating a good fit in each fre- 
quency channel, the merit function is on average approxi- 
mately equal to the number of degrees of freedom, i>, such 
that (x 2 /f) ~ 1-05. Under the assumption that the er- 
rors in the measured Stokes parameters are normally dis- 
tributed, the curvature matrix of equation (8) is inverted 
to yield the covariance matrix of the standard errors in 
each of the model parameters. 

In Figure 3 are plotted the instrumental response pa- 
rameters from one of the CPSR-II bands at the reference 
epoch. The linear dependence between differential phase 
and frequency indicates a signal path length mismatch be- 
tween the two polarizations. The receptor ellipticities have 
an average value of ~ 5.7 degrees, corresponding to a rota- 
tion of the polarization vector about the [/-axis by ~ 0.2 
radians. Consequently, the degree of mixing between lin- 
ear and circular polarizations reaches approximately 20% 
for uncalibrated signals with polarization vectors lying 
near the Q — V plane. Clearly, this level of distortion 
cannot be treated as a second-order effect. 

The mean polarimetric pulse profile of PSR J0437— 4715 
is plotted in Figure 4. Here, the position angle is equal to 
9 and the colatitude is given by 7r/2 — 2e. When com- 
pared with cylindrical coordinates, as plotted in Figure 3 
of van Straten (2002), spherical coordinates offer a num- 
ber of conceptual advantages. For instance, both the total 
and polarized intensities may be plotted on the same log- 
arithmic scale, thereby enhancing the smaller features of 
these profiles. In addition, the chosen normalization al- 
lows the statistical significance of the measured quantities 
to be approximated without explicitly plotting the esti- 
mated error bars. Finally, the correlations between rapid 
transitions in position angle and dips in polarized intensity 
are more obviously apparent, for example, at pulse phases 
approximately equal to 0.26, 0.50, and 0.84. These regions 
are interpreted as transitions between two highly polar- 
ized, nearly orthogonal, superposed modes of radiation, 
possibly the natural modes of the pulsar magnetospheric 
plasma (Pctrova 2001, and references therein). 



6. CONCLUSIONS 

A fundamentally different approach to modeling the re- 
ception and propagation of polarized radiation has been 
presented in this paper. Based on the formalism devel- 
oped in Sections 3 and 4, a detailed calibration model has 
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been constructed from simple, modular components, in- 
corporating multiple signal paths and source polarizations 
as constraints. Many of the limiting assumptions made in 
previous treatments have been eliminated in this analysis, 
enabling its application in a wider variety of experiments. 
This development is increasingly relevant in the context 
of design considerations for the next generation of radio 
telescopes and instrumentation. For instance, it has been 
proposed that the Square Kilometer Array (SKA) will be 
an interferometric array of a large number of low-cost an- 
tennas. Especially in the case of fixed dipole pairs, the po- 
larizations of the receptors may be highly non-orthogonal; 
therefore, it is important to develop more sophisticated 
theories and techniques of polarimetric calibration. 

Due to the current lack of a standard catalog of multi- 
frequency polarimetric pulse profiles, it is difficult to con- 
firm the validity of any calibration technique. Future ef- 
forts should include the establishment of a set of well- 
calibrated sources, including a number of stable pulsars 
that are monitored on a regular basis at more than one 
observatory and at a variety of radio frequencies. Espe- 
cially at lower frequencies, other propagation effects such 
as those arising in the ionosphere will have to be included 
in the reception model applied to these observations. Once 
established, a relatively short integration on a bright, cal- 
ibrated pulsar could be used to quickly determine the in- 
strumental response without the need for time-intensive 
techniques such as the one described in this paper. 

In order that the results presented in Section 5.2 may 
be used to calibrate other observations, they are stored 
as binary table extensions of a PSRFITS file (Hotan, van 
Straten & Manchester 2004), a pulsar data format defined 
using the Flexible Image Transport System (Hanisch et al. 
2001). Both the PSRFITS file format and its associated vi- 
sualization and reduction software have been openly devel- 
oped in an effort to facilitate the exchange of pulsar astro- 



nomical data between observatories and research groups. 

As the current work is based on the Jones calculus, it is 
limited to describing only "non-depolarizing" or "pure" 
component transformations; it is not possible to repre- 
sent the conversion of a fully polarized signal into a par- 
tially polarized one using Jones operators (Hamaker, Breg- 
man & Sault 1996). This is demonstrated by noting that 
det(AB) = dct(A)dct(B) and, for a fully polarized sig- 
nal, det(p) = 0. Certain components of the instrumental 
apparatus may not be pure and therefore may require a 
Mueller matrix in order to model their effect. For example, 
the process of two-bit quantization may act to depolarize 
strong input signals, an effect that has yet to be studied in 
rigorous detail. Other phenomenon, such as bandwidth de- 
polarization, may not require a Mueller matrix description 
if they can be treated by the application of phase-coherent 
matrix convolution (van Straten 2002). 

The formalism developed in this paper may be equally 
well applied to the polarimetric calibration of a phased 
array. In this case, it is not necessary to describe the com- 
plete instrument in terms of the complex gains, orienta- 
tions, and ellipticities of its individual receptors. Rather, 
it is more useful to consider only the total instrumental 
transformation given by the sum of the Jones matrices 
of the individual antennas. This development will be the 
subject of future work. 

The Parkes Observatory is part of the Australia Tele- 
scope which is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO. 
This research greatly benefited from lectures presented 
by Johan Hamaker from 14 January to 18 March, 2003. 
Thanks also to the Swinburne University of Technology 
Pulsar Group for providing the CPSR-II observations and 
to Ben Stappers for advice on the text. 



APPENDIX 

GEOMETRIC INTERPRETATION 

The 2x2 complex matrices with unit determinant may be parameterized by exp[(/3 rh + i(f>n)'(r]. The unit vectors, rh 
and fi, as well as the dimensionless values, (3 and <j>, have meaningful geometric interpretations in the four-dimensional 
space of the Stokes parameters. Under the congruence transformation of equation (3), the Hermitian matrices, 

Brri(/3) = exp(/3m • er) = cosh/3er + sinh/3ra • er, (Al) 

effect a Lorcntz boost of the Stokes 4- vector along the axis rh by an impact parameter 2(3. Likewise, the unitary matrices, 

R-ft(0) = exp(i<j>fi • er) = cos</>er + ismcfin • er, (A2) 

rotate the Stokes polarization vector about the axis ft by an angle 2<f). In applying this geometric representation, the 
following conventions are used. First, the Pauli spin matrices are defined as 

Second, the three-dimensional space of the Stokes polarization vector, S = (Si, S2, S3), is spanned by the orthonormal 
basis vectors, Sfc, such that Jfc • S = Sk and §f. • a — erk- Finally, Stokes Q, U, and V are calculated by the projection of 
S onto the Stokes unit vectors, q, u, and v, respectively. This relationship is summarized by p — (Q, U, V) = R T S, where 
R = (q u v) is a three-dimensional rotation matrix with columns defined by the Stokes unit vectors. The orientation of 
these basis vectors with respect to Sf. depends upon the reference frame in which the electric field vector is represented. 
In the Cartesian basis, the plane wave propagates along the z-axis, the electric field is decomposed into its projection 
along the x and y axes, and R is equal to the 3x3 identity matrix. 
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Table 1 
Signal Path Transformations 



Source 


Transformation 


Noise Diode 
PSRJ0437-4715 
3C 218 


J(i)=GB Sl ( 7 )R Sl (#)C 

J(t)R«(*) 
G H B Sl ( lH )H Sl ( VH )C 



<&Q 




1.2 
1.15 

1.1 
1.05 



'_ I 1 






- 1 , 


i i 




- 



-100 -50 50 

Parallactic Angle (degrees) 



100 



Fig. 1. — Observed Stokes parameters from PSR J0437— 4715, normalized by the invariant interval and plotted as a function of parallactic 
angle. These data correspond to a pulse phase of approximately 0.485 in Figure 4, as observed in a single 500 kHz channel centered at 
1324.75 MHz. The Stokes parameters predicted by the best-fit model are drawn with solid lines. 




1340 
Frequency (MHz) 

Fig. 2. — Stokes parameters of the noise diode reference signal, plotted as a function of observing frequency. The modeled values of 
Stokes Q, U, and V [C\, C2, and C3) are specified as percentages of the reference flux density, Cq. Error bars correspond to the formal 
standard deviations of the model parameters derived from the best-fit covariance matrix. 
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Fig. 3. — Best-fit model parameters as a function of observing frequency. From top to bottom are plotted the orientation of receptor 1 
with respect to receptor 0, 0\, the cllipticitics of the receptors, e*,, the differential phase, <j> = <p(0), the differential gain, 7, and the absolute 
gain, G, specified in units of the square root of the reference flux density, co = \J Co- In the panel showing the cllipticitics, black and grey 
correspond to receptors and 1, respectively. As in Figure 2, error bars correspond to the formal standard errors of the model parameters. 
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Pulse Phase 

Fig. 4. — Mean polarization of PSR J0437— 4715, plotted as a function of pulse phase using polar coordinates: orientation, 6, cllipticity, e, 
and polarized intensity, S = \S\ (plotted in grey below the total intensity, So). Flux densities are normalized by <jq, the r.m.s. of the off-pulse 
total intensity phase bins. Data were integrated over a 64 MHz band centered at 1341MHz for approximately 8 hours. 
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DEGENERACY UNDER COMMUTATION 

The Pauli spin matrices satisfy <Ji (Tj = 5ij er + iajk er fc , where Cijk is the permutation symbol and summation over 
the index k is implied. Using this identity, the product of two arbitrary matrices, A and B, may be expanded to yield 

AB = (a er + a • er) (b er + b • er) = (ab + a • b)er + (ab + ba + i a X b) -cr. (Bl) 

In this form, it can be seen that A and B commute if and only if a and b are collinear, enabling simple statements to be 
made regarding the uniqueness of solutions to the polarization measurement equation. 

For example, consider the parallactic angle rotation given by R{;( ( f>) = cos$ero + isin${) • er. As illustrated by 
equation (Bl), any matrix of the form, V = Vq er + Vv • er, commutes freely with R«($). Therefore, if both J and p m 
satisfy equation (20), then 3 V = JV 1 and p v m = Vp m V^ are also solutions to this equation. This degeneracy exists 
regardless of the parameterization of J or p m , proving that no unique solution to the polarization measurement equation 
may be derived solely from observations of unknown sources at multiple parallactic angles. 

The consequences of this degeneracy are illustrated by noting that V is unimodular and may be decomposed into a 
boost along the v axis, B«(/3„), and a rotation about this axis, Re ((/>„). That is, unless V is determined through some 
other means, there remains an unknown position angle error as well as an unknown degree of mixing between the total 
intensity and circular polarization. An observation of a (possibly unpolarized) source with a well-determined degree of 
circular polarization may be used to determine f3 v ; additionally, an observation of a linearly polarized source with a known 
position angle (such as the noise diode included in many receiver packages) may be used to constrain <ft v . 
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